The mechanism of ion induced amorphization in Si 
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Some results on damage build up in, and amorphization of, Si, induced by 25-30 keV Al^, Si^ 
and Cs~ ions, at room temperature, are reported. We show that at low energy, amorphization is a 
nucleation and growth process, based on the direct impact mechanism. With an Avrami exponent 
~ 1.6, the growth towards amorphization seems to be diffusion limited. A transition to a completely 
amorphized state is indicated at a dose exceeding 17 eV/atom, which is higher than 6-12 eV/atom 
as predicted by simulations. The observed higher threshold could be due to temperature effects 
although an underestimation of keV-energy recoils, in simulation, may not be ruled out. 



The study of ion implantation induced damage and 
recovery is a very important area of research in semi- 
conductor processing, particularly involving Si. In all 
the cases involving doping through ion implantation, a 
damage layer is formed which must get back to a defect 
free crystalline state for any later application. In view of 
this, defect production by ions, together with its growth 
and annealing behavior, constitute an important area of 
study. 

Ion implantation induced amorphization has been the 
subject of intense research in which there is a long stand- 
ing debate dating back to the seventies [l|, . One of the 
view points is that amorphization is caused from over- 
lapping of amorphized pockets formed from defects cre- 
ated by individual ion cascades or direct impact. This is 
the so called heterogeneous amorphization as suggested 
by Morehead and Crowder [31]. There are experimen- 
tal data in support of this [1, H, 0]. Competing with 
this, there is a homogeneous amorphization mechanism 
where passage of the energetic ions results in the forma- 
tion of a large number of bond defects (consisting of iso- 
lated point defects and interstitial- vacancy complexes). 
During implantation, with formation energies of ~3 eV, 
these defects can be uniformly produced in the system 
and when their concentration increases beyond a certain 
limit the lattice becomes unstable resulting in a collapse 
to an amorphized state ■ There are experimental data 
in support of this as well [1, H, [13, [HI which also sug- 
gest amorphisation in Si is more like a phase transition 
[12j induced by an accumulation of a sufficient number 
of defects. Reference ^1] provides a recent review on the 
subject. 

In this letter we present some results regarding dam- 
age production and growth from low energy Al^, Si^ 
and a similar mass Cs^ ions in Si at room temperature. 
Transmission electron microscopy (TEM) and channeling 
Rutherford backscattering spectrometry (RBS/C) have 
been used for sample characterization. Increase in flu- 
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ence, 0, results in an increase in the number of atoms 
displaced from lattice positions leading to an increase in 
backscattered yield in the aligned condition. The growth 
in amorphous fraction, fa, is obtained through a relative 
growth in the surface peak intensity, over and above that 
of a virgin sample. Three different growth models, all 
based on the direct impact mechanism, have been used 
to explain the data. We show, at low energy, amorphisa- 
tion proceeds via a nucleation and growth process under 
the direct impact mechanism. The growth against dose 
is consistent with a Kolmogorov-Johnson-Mehl- Avrami 
(KJMA) equation with an exponent ~ 1.6 as against 
a value between 3 and 4, predicted for 3D growth with 
a constant or slowly varying nucleation rate [13i] . This 
suggests the growth to be diffusion limited. Complete 
amorphization is found to occur above a dose, Z?, of 17 
eV/atom, much higher than 6-12 eV/atom as suggested 
by Molecular Dynamics (MD) simulations [2]. At doses 
below the saturation level, the amorphous to crystalline 
(a/c) interface is found to be rough. Since the experi- 
ments were carried out at room temperature, the rough 
interface can not be due to immobile defects as observed 
in case of a similar study [l3| carried out at liquid N2 
temperature. However, at ~ 50 eV/atom a completely 
relaxed amorphous phase, with a smooth a/c interface, 
has been observed. 

The sample preparation involved both cluster (AI5 
and Si5) and single atom (Cs) implantation into Si(lOO) 
(p-type, 20 rJcm) substrates, at room temperature, at 
low beam currents of ~ 2-3 nA. Five samples were im- 
planted with Sig" ions, at 5 keV/atom, to cluster fluence 
of 2 X 10", 4 X 10", 4 X 10^2, 1.2 x lO^^ and 1 x lO^'' 
cm^^ respectively. Two more samples were implanted 
with 25 keV Cs~ ions and 30 keV Al^ ions (6 keV/atom) 
to fluence of 6 x 10^'^ cm~^ and 1.6 x 10^^ cm~^ respec- 
tively. Excepting few selected ones, subjected to high 
resolution (HR) cross sectional TEM (XTEM) imaging 
(at 200 keV), all the samples were subjected to RBS/C 
measurements using 1.35 MeV He'^ ions at a scattering 
angle of 130°. All the implantation and measurements 
were carried out using the facilities at the Institute of 
Physics (lOP), Bhubaneswar. For quantitative analysis, 
the amorphized fraction, fa, for a given ffuence 0, was de- 
termined from the surface peak intensity in RBS/C data. 
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FIG. 1: Growth of amorphous fraction, /a, with implantation 
dose, D, in eV/atom. The fitted functions correspond to two 
of the models as discussed in the text. KJMA fits to the 
present data and data extracted from simulation results of 
Ref [2], are shown in the inset. The three fits yield Avrami 
exponents, n, close to 1.6 as indicated. A schematic diagram 
of the RBS/C experimental geometry is shown in the top right 
corner. 



relative to a virgin sample. This value was normalized to 
unity at the highest fluence (corresponding to complete 
amorphization). Since AI5, Sis and Cs are different sys- 
tems we convert the fluence, to dose, D (eV/atom). 
This is done using the equation, D — cf)E/{RpN), Rp, E 
and N representing the projected range of the implanted 
atom, its energy and the atomic density of the matrix 
(5 X 10^^ cm~^) respectively. Since the clusters break 
almost immediately upon impact, the damage produced 
due to their implantation is mainly due to 5 or 6 keV Si or 
Al atoms respectively which have ranges of 10.2 and 12.7 
nm in Si [l3|. Based on TEM data (shown later), for Si 
and Cs atoms, we take the ranges to be 10 nm and 20 nm 
respectively. The results are shown in Fig. 1. The 95% 
level in fa is seen to occur at a dose of ~ 17.5 eV/atom, 
complete saturation occurring around 20 eV/atom. This 
is higher than the saturation value of 6-12 eV/atom as 
indicated by simulations 2]. 

The results can be explained in terms of a combined 
Direct impact and defect stimulated mechanism [l^ . This 
is mainly because presence of defects is known to result 
in a faster growth in the amorphized fraction jT^, ITl| . 
The cross-section for the defect stimulated process, as, 
depends on the same for the direct impact process, Ga- 
In such a case, fa = l-{aa+crs) I {crs+craexp{{aa+cFs)D)) 
[Hi. A fit to the data yields values of 0.254 (± 0.012) 
and 0.04 (± 0.005) for ct^ and da respectively. To 
look at cascade overlap effects we have used the Cas- 
cade quenching and recrystallization model proposed by 
Wang et al [17|. Here amorphization is assumed to 
be the result of direct impact of energetic ions, how- 
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FIG. 2: (a) XTEM images of the Si sample irradiated with Sis 
to a cluster fluence, 0, of 4 x 10^^ cm~^; (b) an HR image of 
the undamaged region as marked in (a); (c)&(d) show two HR 
images of two marked boxed in the defected region, showing 
(within circles) an amorphous patch and dislocations, within 
the ion EOR. 




FIG. 3: (a) XTEM images of the Si sample irradiated with 
Sis clusters to a fluence, 0, of 1.2 x 10^^ cm~^; (b) same as (a) 
at a lower resolution; (c) a Fourier transform of the region as 
marked by the box in (a); (d) the filtered image of the marked 
region. The circles indicate representative amorphous patches 
present in the damaged region. 



ever, under a competing effect of partial recrystalliza- 
tion. Using a recrystallization efficiency. A, one can 
write fa = I - I/^JA (1 - A)exp[2{l - A)Dn], where 
Dn = kD, represents a dimensionless normalized dose. 
The parameter k is related to the cross sectional area of 
the damaged region at any given dose. A fit using the 
above model yields a value of 0.893 (±0.002) for A, to- 
gether with a k value of 2.214 (±0.048). This k value 
indicates, on an average, a 5 keV Si atom generates dam- 
age in a cylindrical region of cross sectional area of 22 
nm^ p^ . corresponding to a track radius of 2.65 nm. 
Now we look at our fit values against the simulation re- 
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suits of Nord et al 0] (as given in their table-IV). In 
our case saturation occurs at 20 eV/atom indicating a 
normalized dose, D„ of about 44. With almost identi- 
cal A values, simulations for the two cases viz with 3 
eV - 2.0 keV and 15 eV - 2 keV recoils, show k values 
of 3.822 and 6.326 respectively. The corresponding D 
values of 12 eV/atom and 6.5 eV/atom result in very 
similar Dn values of 45.86 and 41, in good agreement 
with present data. A _D„ value of 44 indicates a com- 
plete amorphisation at a (/> of 4 x 10^* atoms-cm~^ in 
good agreement with earlier data [l3|- The growth in 
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FIG. 4: (a) An XTEM image of the Si sample irradiated with 
Cs atoms to a fluence, 0, of 6 x lO'^^ cm~^; (b) same as (a) 
at a lower resolution; (c) a Fourier transform of the region as 
marked by the box in (a). The circles show some representa- 
tive crystalline patches in a largely amorphized medium. 
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FIG. 5: (a) An XTEM image of the Si sample irradiated with 
Sis clusters to a fluence, 4>, oi 1 x 10^* cm~^; (b) same as (a) 
at a lower resolution; (c) a Fourier transform of the region as 
marked by the box in (a). 



fa with dose, D, can also be studied using a Nucleation 
and growth model based on the KJMA equation ,13] of 
the form fa = ^ — exp[—RG"~^D"^]. Here R is related to 



the nucleation rate (dependent on a a) while G is related 
to the growth velocity (dependent on ct^), at the interface 
[l8l |. The Avrami exponent, n, as determined from the 
slope of ln\—ln{\ — fa)] vs InD (Fig. 1) yields a value 1.6 
±0.06 (close to 3/2), indicating a diffusion hmited nucle- 
ation and growth mechanism to be valid (2^ 1 . Simulation 
results of Nord et al, [3| , for realistic (variable energy) re- 
coil distributions, yield very similar n values (Fig. 1) in 
agreement with our findings. This is in contrast to some 
earlier results on Ge induced amorphisation of Si carried 
out at higher energies, where n was shown to be 3.5 [2l| . 
Since n has almost the same value in both simulated and 
experimental data, different as values indicates unequal 
growth velocities at the a/c interface. 

Fig. 2 (a) shows an XTEM picture of the defected 
region as obtained for an Sis fluence of 4 x lO'^^ cm~^, 
corresponding to a dose, D, of 0.2 eV/atom. Fig. 2(b) 
corresponds to an undamaged region. Fig. 2(c) & (d) 
show defects in the form of an amorphized patch and dis- 
locations respectively, present within a depth of 10 nm 
(range of 5 keV Si in Si) . HR XTEM images for the sam- 
ple implanted with Sis fluence of 1.2 x 10^^ cm~^ (dose 
~ 6 eV/atom) are shown in Fig. 3 (a) & (b). In this case 
a heavily damaged region, extending from the ion end- 
of-range (FOR) up to the surface, can be seen. One can 
notice a rough interface together with a lot of crystalline 
patches in the near surface region. This is due to spatial 
overlap of collision cascades at the ion EOR where a high 
density of low energy recoils is expected to be produced. 
A Fourier transformed (FT) image of a 10 nm size region 
(Fig. 3 (c)) shows diffused < 111 >-diffraction spots 
along with shifted lobes indicating incomplete amorphi- 
sation and strain in the system. A Fourier inversion of 
the diffraction image (Fig. 3(d)) clearly shows the heav- 
ily damaged lattice with misoriented planes. One can 
also see some amorphized patches (shown in circles) in 
the damaged region. Compared to this, the XTEM im- 
ages for the Cs implantation case (Fig. 4 (a) & (b)), for 
a dose of ^ 15 eV/atom, show a greater degree of amor- 
phization. The FT image (Fig. 4 (c)), of a damaged 
patch as marked, shows a continuous ring indicating a 
largely amorphized layer. The direct image also shows 
this, however, with some crystalline patches still present. 
The corresponding fluence for 5 keV Si in Si, 1.5 x 10^"* 
cm"^, is very close to the threshold for amorphisation 
(the 95% level occurs at 17.5 eV/atom). But the inter- 
face is still very rough. HR XTEM images of the dam- 
aged region obtained with an Sis implantation fluence of 
1 X 10"'^'' cm~^ are shown in Fig. 5 (a)&(b). This case 
corresponds to a dose of about 50 eV/atom which is well 
above the amorphization threshold. Fig. 5(a) shows com- 
plete amorphization, the diffraction image (Fig. 5(c)) of 
a representative marked region showing two clear rings. 
The a/c interface is also seen to be rather flat, mainly 
due to stress relaxation [2^. 

The above observations can be explained as follows. Si 
(or Al) atoms at 5 (or 6 keV) or 25 keV Cs atoms, in 
addition to creating point defects, produce amorphous 
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patches in the implanted region. Some of these, are 
seen in Fig. 2(c). With increase in dose, with larger 
overlap of collision cascade regions there is a growth in 
the amorphized volume extending from the EOR towards 
the surface. Fig. 3 and Fig. 4 show this, with chunks 
of crystalline material present near the surface. These 
crystalline regions finally get amorphized from impact of 
low energy recoils due to successive ion passes. However, 
the growth in amorphisation is diffusion limited as the 
Avrami exponent, n, has a value close to 3/2 [l^. We 
believe, the above diffusion limited growth mechanism 
also holds for amorphisation near the EOR. In this case 
the rough interface is smoothed out because of the depth 
limit imposed by the ion energy in addition to stress re- 
laxation. 

For most of the potentials used for simulation of the Si 
structure, an energy deposition of 6-8 eV/atom is needed 
for getting an amorphous structure This is small 
compared to the amorphizing dose (> 17 eV/atom), as 
obtained here. One of the reasons could be that the ex- 
periments were carried out at room temperature {300K) 
where small amorphous patches can get recrystallized 
even without any ion passage [25j . leading to a smaller 
growth rate. Another reason could be a possible under- 
estimation of the recrystallization effects due to keV en- 
ergy recoils in the simulation. At higher beam energies 



and higher beam currents, amorphisation will be sup- 
pressed by recrystallization induced by relatively larger 
number of recoils, pushing the amorphization threshold 
to a higher dose. This could be the reason behind get- 
ting an amorphization threshold of 4 x 10^'' cm~^, with 
230 keV Si IQ], as compared to almost half that in the 
present case. Beam energy and current induced anneal- 
ing effects may even result in recrystallization of amor- 
phous patches leading to the formation of point defects 
and related complexes [l^, ■ In case the concentration 
such defects becomes too high homogeneous amorphiza- 
tion may take place. This does not happen in the present 
case. 

To conclude, using 30 keV AI5, 25 keV Si^ and 25 keV 
Cs^ ions, we have shown that ion induced amorphisation 
in Si, at low energy, is very much dependent on the direct 
impact mechanism resulting in heterogeneous amorphi- 
sation. At these energies, complete amorphisation seems 
to take place through a diffusion limited, nucleation and 
growth process. 
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